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Abstract 

We study the dynamics of a one-dimensional non-linear and non- 
local drift-diffusion equation set in the half-line, with the coupling 
involving the trace value on the boundary. The initial mass M of the 
density determines the behaviour of the equation: attraction to self 
similar profile, to a steady state of finite time blow up for supercritical 
mass. Using the logarithmic Sobolev and the HWI inequalities we ob- 
tain a rate of convergence for the cases subcritical and critical mass. 
Moreover, we prove a comparison principle on the equation obtained 
after space integration. This concentration-comparison principle al- 
lows proving blow-up of solutions for large initial data without any 
monotonicity assumption on the initial data. 

Keywords: Cell polarisation, global existence, blow-up, asymptotic 
convergence, entropy method, Keller-Segel system, logarithmic Sobolev 
inequality, HWI inequality. 

1 Introduction 

In this paper we improve the analysis of a one-dimensional non-linear and 
non-local convection-diffusion equation introduced in a previous paper [3] 

dtn{t, x) — dxx'n{t, x) = n{t, 0)dxn{t, x), x > 0, 

dxnit,0)+n{t,Of = 0, (1) 

n{0,x) = n^{x), X > 0. 

The boundary condition ensures mass conservation: 

n{t,x)dx= / n°(x)(ix = M. (2) 

^0 
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The previous model ([T|) comes from the models given in [9] and [3] under 
a radial symmetry assumption. These latter models describe cell polarisa- 
tion. Cell polarisation is a major step involved in several important cellular 
processes such as directional migration, growth, oriented secretion, cell divi- 
sion, mating or morphogenesis. When a cell is not polarised proteins Cdc42 
are uniformly distributed on the membrane while polarisation is character- 
ized by a concentration of proteins in a small area of the cell membrane. 

In [9] and [3j, in dimension higher than two, a class of non-linear convection- 
diffusion models were designed, and studied, for cell polarisation. In these 
models there is a coupling between the evolution of proteins and the dynam- 
ics of the cytoskeleton: the proteins diffuse and they are actively transported 
along tubes or filaments towards the membrane. The advection field is ob- 
tained through a coupling with the concentration of markers on the mem- 
brane. The resulting motion is a biased diffusion regulated by the markers 
themselves. 

Of special interest is the fact that solutions of ([T]) may become un- 
bounded in finite time (so-called blow-up). Blow-up of solution of ([T]) means 
that convection wins over diffusion. In such a situation markers concentrate 
on the membrane of the cell. We recall this result from |3]. 

Theorem 1.1. Assume M > 1. Any weak solution of equation ([TJ (in the 
sense of Definition \2. 1\ ) with non-increasing initial data blows-up infinite 
time. 

In this paper, we state a so-called concentration-comparison principle on 
the equation obtained after space integration. This principle together with 
the use of a self-similar supersolution gives estimates from above on n{t, 0) 
and a comparison with a suitable heat equation allow obtaining estimates 
from below on n(t, 0) and extending the blow-up result to any initial data 
above the critical mass. 

Proposition 1.2 (Concentration-comparison principle). Letni,n2 G C^{0,T;M.-^-) 
be two solutions of equation ([T]) satisfying 



On the other hand, global existence and asymptotic behaviour in the 
sub-critical and critical cases, M < 1, were established in [3]. 

Theorem 1.3. Assume that the initial datum satisfies both G 
^-"^((1 + x)dx) and 




Vx > 0, 



then, we have 
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n^{x) (log 'nP{x))+dx < +00. Assume in addition that M < 1, then there 
exists a global weak solution that satisfies the following estimates for all 
T>0, 



sup 

te(o 

T roo 



roc 

p / n{t, x) {log n{t,x))- 
,T) Jo 



dx < +00 , 



Jo 



n{t,x) {dx log n{t,x))'^ dxdt < +00. 



In the sub-critical case M < 1 the solution strongly converges in towards 
the self-similar profile Ga given by \1S\) in the following sense: 



lim 



n{t, x) 



1 



--G, 



VT+2i \VTT2i 



In the critical case M = 1, assuming in addition that the second momen- 
tum is finite x'^n^{x)dx < +00, the solution strongly converges in 
towards the unique stationary state i'a{x) '■= aexp(— ax), where = 
xnP{x)dx. 

In [3], in the sub-critical case, the asymptotic result was obtained through 
the convergence to zero of a suitable Lyapounov functional L defined by ([13]). 
Here, using a logarithmic Sobolev inequality with a suitable function, we ob- 
tain an exponential decay to equilibrium in self similar variables replacing 
in particular the former result by 

Proposition 1.4. Under the assumptions of theorem we have 



n{t, x) 



1 



zGa 



VTT2i \VT+2i 



< 



c 



(l + 2t)3/2' 



Finally, in the critical case M = 1, we actually improve the result given in 
Theorem 11.31 by precising the speed of convergence. Firstly, we give a rate. 
Secondly, using the HWI inequality, we improve the rate of convergence 
when the third momentum is initially finite. 

We end this introductory Section with one open question that we are 
not able to resolve: obtain a blow-up profile for large initial datum. 

The plan of this work is the following. First, we state a concentration- 
comparison principle and the finite time blow-up for supercritical mass. 
Then, we give a quantitative argument to the systematic blow-up for su- 
percritical mass. In a third step we study the sub-critical mass and we give 
a rate to the self similar decay. Finally, we study the critical mass, by using 
first a Lyapounov approach and then the HWI inequality. 
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2 A brief account of some useful facts 

We first recall some facts concerning ([1]), see [3J for more details. Let n be 
a classical solution of ([1]) on (0,T), a straightforward computation of the 
evolution of the entropy yields 



— / n{t,x) log n{t^x)dx = — {dxn{t,x) + n{t,0)n{t,x)) 
dt Jo Jo 



dxn{t,x) 
n{t, x) 

n(t, x) {dx log n(t, x)f dx + n{t, Of (3) 
Moreover a proper definition of weak solutions, adapted to our context 



IS 



Definition 2.1. IS] We say that n{t,x) is a weak solution of ([T]) on (0,T) 
if it satisfies: 

n G L°°(0,T;L^(M+)), d^n € L^((0,r) x M+) , 
and n{t,x) is a solution of ([T]) in the sense of distributions in 2?'(M-|_). 

Since the flux dxn{t,x) + n{t,0)n{t,x) belongs to L^((0, T) x M+), the 
solution is well-defined in the distributional sense under the assumptions of 
definition (|2.ip . In fact we can write n{t, 0)dt = — dxn{t, x)dxdt. 

Prom the maximum principle [7] it follows that n > if n*^ > 0. Fur- 
thermore weak solutions in the sense of Definition 12.11 are mass-preserving: 

M= nP{x)dx = / n{t,x)dx. 
Jo Jo 

Throughout the paper, for n > G L^(M+), we define the generalized 
distribution function in the following way: 

N{x) = / n{y)dy, 
Jo 

and for any weak solution n in the sense of definition 12. H the first momentum 
is defined as 

poo 

J{t) = / xn{t, x)dx. 



Let u{x)dx and u{x)dx be measures with smooth densities on M_(., we note 
the relative entropy 



u{x) log ( ) dx 



u[x 
u{x) 

and the Fisher information 

2 



I{u\u)= I u{x){dx{log^] ] dx. 
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We also define the quadratic Wasserstein distance between two probability 
measures /i and v on with finite second momentum as in [16], 



W2{li,v)= inf . |x-y|2(i7r(x,y). (4) 

where n(/i, v) denotes the set of probability measures on x with 
marginals ^ and i.e. such that for all test functions (p and m. a, suitable 
class of test functions, 

{4>{x) +i){y))dTT{x,y) = 4>{x)dfi{x) + tp{y)dv{y). 

xIR+ JO Jo 

In particular, if is a Dirac ineasurG 8a with G M-i- we havG 

/"OO 

W2{^i,5af= / |x-a|2(i/i(x). (5) 

JO 

A simple computation on the first momentum (if it is initially finite) leads 
to 

|j(t) = n(t,0)(l-M). (6) 
We recall the three situations: 

• when M = 1, there is a balance between drift and diffusion and we 
expect a convergence to a steady state (see section [5]), 

• when M < 1, the diffusion drives the equation and we have a self 
similar behaviour see section HI 

• when M > 1 , solution blows up in finite time see section [31 

Remark 1. Such a critical mass phenomenon (global existence versus blow- 
up) has been widely studied for the Keller-Segel system (also known as the 
Smoluchowski-Poisson system) in two dimensions of space, see [7^ e.g. 
and the references therhein. The equation represents in some sense a car- 
icatural version of the classical Keller-Segel system in the half-line (0, +oo). 

Remark 2. There is a strong connection between the equation under interest 
here ([ip and the one- dimensional Stefan problem. The later writes 1 11^ : 

dtu{t, x) = dxxu{t, x) , t> 0, X e (-00, s{t)) , 

lim^^_oo d::,u{t, X) =0, U{t, s{t)) = , d::,u{t, s{t)) = -s'{t) . 

The temperature is initially non-negative: u{0,x) = uo{x) > 0. By perform- 
ing the following change of variables: (f>{t.,x) = —u{t,s{t) — x), we get an 
equation that is linked to (QP by n{t, x) = dx4>{t, x). This connection provides 
some insights concerning the possible continuation of solutions after blow- 
up fll]/ . This question has raised a lot of interest in the past recent years 
110[ {Bjj. It is postulated in fi i]/ that the one- dimensional Stefan problem is 
generically non continuable after the blow-up time. 
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3 Concentration-comparison principle and finite 
time blow up for supercritical mass 

The main result of this section is a concentration-comparison principle given 
by Proposition 11.21 that we first prove. Then, we use this principle to com- 
pare solutions of equation ([1]), rewritten as an equation on A^, to both sub- 
and supersolutions. Finally we provide a quantitative blow-up argument by 
estimating the blow-up time T* defined by J(T*) = 0. 

3.1 Proof of Proposition [TT2] 

We actually prove the following stronger lemma: 

Lemma 3.1. Let N, N,N_ be nondecreasing (in space) functions in C^(0, T; C^(R-i-)) 
satisfying 

d^N,d,N,d^NeL^{0,T;L^{R+)) 

and 

' dtN{t,x) - d^^N{t,x) - d.,N{t,0)d,N{t,x) = 0, 
^ dtN{t, x) - d,,N{t, x) - d,N{t, 0)d^N{t, x) > 0, 

dtN{t, x) - d^^N{t, x) - d.^N{t, 0)d^N{t, x) < 0, 
^ N{t, 0) = N{t, 0) = N{t, 0) = 0. 

Assume that N{0, .) < iV(0, .) < iV(0, .), then for < t < T we have that 

N{t,-)<N{t,.)<N{t,.). 

Proof. It is obviously sufficient to prove the comparison N > N_. Firstly, we 
assume that 

d^N{0,0) > 9,iV(0,0) (> 0). 

Integrating equation JT]) in space, we see that 5N = N — N_ satisfies the 
parabolic equation 

dt6N{t, x) - d.JNit, x) - d,N{t, 0)d,6Nit, x) = (9,iV(t, 0) - a,iV(t, 0))9,iV(t, x), 
Ni{t,0) = 0. 

Since we supposed dxN{0,0) > dxN,{0,0), it remains true at least until 
a time T' s]0, T[, we choose the maximal T' possible. On the time interval 
[0, T'] we have, since dxN. > 0, 

dt6N{t,x) - dxx6N{t,x) - dxN{t,0)dx5N{t,x) > 0, 
Ni{t,0)=0. 
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Hence > iV on ]0, r'[xM^ by strong maximum principle [7J. Furthermore, 
by Hopf Lemma (see [7j ) , we also have 

dJNiT',0) = d,N{T',0) - d,N{T',0) > 0. 

As T' is maximal we immediately conclude that T' = T. To treat the case 
da:N{0,0) = d^N{0,0), we use the following 

Lemma 3.2. Let N be a supersolution as in lemma \3.1l denote K = 
sup[o,r] , then 

N\t, x) = N{t, x) + ee^^'\l - e"^-^^) 
is a supersolution on [0,T] in the sense of lemma CO] for < e < ^e~^^^'^ . 

Proof. It is straightforward that dxN'^ > and -/V'^(i,0) = 0. We essentially 
have to prove the inequation 

p = dtN^ - d^^N^ - d^N^it, 0)dxN^ > 0. 

We denote 

N'{t,x) = N{t,x) + P{t){l - e-°^). 

Conditions of the lemma ensure P{t) < i for < t < T. Since is a 
supersolution, we have 

r > /3'(t)(l-e-°^)+a2/3(t)e-"^-/3(t)a9^.A^-/3(t)ae-"^5^.iV(t, 0)-a2/32g-ax_ 
This leads to, by definition of K, 

r>(3'- (3aK + /3e-"^ j + a'^ {1 - (3) - aK^ 

and the right hand side is nonnegative. We choose the = ee"^*, so that 
r > /3e-°^ (a2(l - ee"-^^) - 2aK) , 

and finally, choosing a = AK and e < - — ^ — we have /"^ > 0. □ 
End of the proof of lemma \3.1\ 

Thanks to the first step of the proof, since for e > we have dxN'^{t, 0) > 
dxN{t, 0) > dxN{t, 0), we can compare N'{t, x) > N{t, x) for < t < T and 
e small enough. Letting e — )■ we can conclude N{t,x) > N.{t,x). □ 
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3.2 Estimations on n{t, 0) and local existence for supercritical 
mass 

Classically, the concentration-comparison principle allows us comparing so- 
lutions of d?]) to sub- and supersolutions. Since the proof is very much alike 
the one for the comparison principle (also based on Hopf Lemma) we do not 
repeat it. 

The main use of this Lemma is the following: 

Lemma 3.3 (Self-similar supersolutions). Let A G M and f € C'^(M-|-,M+) 
satisfying 

'f'{x)>0, Vx>0, 
< /(O) = 0, 

, Xxf'ix) - fix) - /'(0)/'(x) > 0, Vx > 0. 
Then the function Nj^x defined by 

is a supersolution of ([7]) for all t > if X < and for t < if X > 0. 

Proof. For 1 — 2Xt > we first notice that 

d 1 _ A 
Vr^2At ~ Vl - 2At^ ■ 

Therefore, by construction, as long as 1 — 2Xt > 0, we see that 

(1 - 2Xt) {dtNxjit, x) - d^^Nxjit, x) - d,Nxj{t, 0)d,Nxj{t, x)) = 

A , f ( , ^ ^ - f" ( , ^ ^ - /'(O)/' ( , ^ ^ > 0. 

VT^^ \y/l-2Xt) \y/l-2Xt) ^ ' Wl-2XtJ - 

Finally, it is easy to prove that such a / exists and that, for all x > 0, it 
verifies ^ 

fix) < /'(O) £ exp (^a|- - /'(0)y) dy. (8) 

□ 
□ 

Remark 3. For X < we can have equality in In such a case, in the 
critical case M = 1, we recover the steady state and in the sub-critical case, 
we recover the self similar profile, see For A > 0, the equality yields 
that /(oo) = +00 which does not correspond to a derivative. Therefore 
a self-similar blow-up profile doesn't seem to appear. 
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This result is particularly useful because it allows comparing a large class 
of initial data to such a supersolution. Indeed, for a good choice of /, by 
comparison principle, we would be able to prove that N{t,x) < Nf^xit,x) 
for t < This will provide an a priori bound on n{t, 0), that is 

Recalling now the entropy 'dissipation' ([3]), for t < and for n'^{x) log n^{x)dx < 
+00, we obtain the following a priori estimates 



sup / n{t, x) {log n{t, x))-^-dx + / / n{t,x){dx log n{t,x)) dx<C{t), 
{o,t) Jo Jo Jo 

which is enough to get compactness, insuring existence of weak solutions 
build as in [4J . Details on existence will be carried out in [13] . To guarantee 
such a bound, we need to be able to compare, that is we need existence of 
X, fj, > such that 



Vx > 



0, ^ /i exp (^a|- - fiy^ dy - N{0, x) > 0. 



3.3 Systematic blow up for supercritical mass : A quantita- 
tive argument 

In [3], the proof of the blow-up in supercritical case is based on the non- 
existence of the first momentum. For the convenience of the reader, we recall 
it now. Assume that M > 1 and that nP is non-increasing, then n(t, •) is also 
non-increasing for any time it exists due to the maximum principle. In fact 
v{t,x) = dxn{t,x) satisfies a parabolic type equation without any source 
term, it is initially non-positive, and it is non-positive on the boundary. 
Therefore —dxn{t,x)/n{t,0) is a probability density at any time t > 0. 
From Jensen's inequality, we deduce that 



n{t,0) J - Jo n{t,0) 

which rewrites as < 2n{t, 0)J{t). Plugging this latter inequality into the 
evolution of the first momentum ([6]) yields that 

J t < J O + ^ / -—ds. 

Introducing next the auxiliary function 

K(t) = J(0) + (1 - M)m2 / 3{sy^ds 

Jo 
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which is positive and satisfies 



d ^^^^ _ (1 - M)M2 1 ^ (1 - M)M2 1 



dt ' ' 2 J(t) - 2 K(t)' 

we deduce that ^ [K(t)^] < (1 — M)M'^, hence a contradiction. 

In this work, we also use the non-existence of decreasing first momentum 
beyond the blow-up time T*, given by J{T*) = 0, but the integrated solution 
gives us a way to find this T* by using equation ([6]): 

d^N{t,0)dt= / n(t,0)dt- ^ ^ 



M -1 



We are interested in using a subsolution of the parabolic equation on N 
to find a blow-up time. A first lower bound on is given by Chebyshev 
inequality. 

Lemma 3.4 (Chebyshev). For allt > and for all x satisfying x > J{t)/M , 
the following inequality holds true 

N(t,x) > M-^. 

X 

Proof. Chebychev inequality applied to the probability distribution n/M 
yields that 

< yn{t,y)dy. 

n 
□ 

We next use this lower bound to define a subsolution of equation ([7]) 
as the solution of a particular heat equation on a bounded domain [0, L]. 
Firstly, recalling that dxN and n{t, 0) are non negative, we deduce that the 
solution N of the classical heat equation on 

'dtN{t,x)-dxxN{t,x) = 0, inM;, 
iV(t,0)=0, (9) 
_iVO = iVO, inM+, 

is a subsolution of equation ([7]). Next, using Lemma ()3.4p . we will consider 
a bounded domain. The advection to disappears in favor of a boundary 
term in x = L. This latter term depends on the decreasing first momentum 
and this gives a weak way to transport to 0. 
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Proposition 3.5. For n solution of equation ([T]) with M >1, for L > 



we have 



n(t,0) > ^-^ 



1+2 5^ i-ir 



exp 



ne N 



n TT 

T 



Furthemore, fort > ^4^) we have 

n{t,0) >k>0. 

Proof. First, if / solution of the bounded heat equation with Dirichlet con 
ditions, 

' dtf{t,x) - dxxf{t,x) = 0, in(0,L), 
< /(t,0) = /(t,L) = 0, 
J%x) = N%x), in(0,L), 
then we recall that 

f{i,x) = ^ sin 



(10) 



nGN* 

with the Fourier coefficients 



nn X 
L 



exp 



n TT 



— ) t)/.. 



Let L be such that L > ^[t^. Since the first momentum J is decreasing, 

from Lemma (fOI) . we deduce that N{t,L) > M - ^ for all t > 0. We 
now define the following problem 



dtN{t,x) -dxxN{t,x) = 0, in(0,L), 
iV(t,0)=0, N{t,L) = M-^ 
^N°{x) = N°{x), in(0,L). 



L ' 



(11) 



We extend the solution N on [L, +oo] by 



N*{t,x) 



N{t,x), ifx<L, 



M 



m 

L ' 



if X > L. 



From Lemma 13.41 together with the fact that the solution of the heat 
equation on M_|_ is a subsolution of equation ([7|) we deduce that A^* is a 
subsolution of equation ([7]). Using Fourier series to solve pi|) . we obtain 
that 



~, , xN(t,L) 
N{t,x) = ^+ ^ sin 

neN* 



nvrx 
L 



exp 



n vr 
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with the Fourier coefficients defined by 

By strong maximum principle |7j, a solution of pop . with initial datum 
A^" > 0, is positive, hence 

~, , xN(t,L) v-^ fmTx\ ( (n'K\'i \ 2N(t,L) , 
iV(t, .) > + 5: sm (^) exp [- (-) t j —1^ (-1)^ 

Furthermore, recalhng that A^(t, 0) = iV(t, 0) = and that iV < on 
[0, +cxd[x [0, L], we deduce that 

n(t,o) = a,iv(t,o) > a,.iv(t,o) > ( i + 2 E (-1)" (" (x)'* 

V n.eN* ^ 

hence, for t > (^)^ log(2), it follows that 

«(«,0)>%a(l-2exp(-g)%))=M<)>0. 

n 
□ 

The previous lower bound on n{t, 0) provides us an upper bound on the 
blow-up time T*. 

Corollary 3.6. For n solution of equation ([1]) with M > 1, we have the 
following upper bound on blow-up time 

^* . -^(0)^ / M-l 
T* <A —TTTT r 1 + 



M2(M- 1) V 6 
Proof. From proposition 13.51 it follows that 



M-l 



ngN* 

(-1)" 



Moreover, recalling that = —j2 , ^or all L > we see that 

- M -Im 6 ■ 
For L optimal, i.e. given by L = ^^7^ — obtain that 

- M2(M-r ' 
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□ 
□ 

Remark 4. The upper bound on the blow-up time T* given in the previous 



this work, convection is described by the decreasing of the first momentum, 
then we have only used diffusion phenomenon in equation (|lip . This could 
explain the difference in results. 

4 Subcritical mass and self similar decay 

For the case M < 1, global existence has been proved in [3]. Here, we are 
interested in the asymptotic behaviour and self similar decay has also been 
exhibited in [3]. If we perform the following change of variable: 



then, the density u{t, y) satisfies 

dru{T, y) - dyyu{T, y) - dy{yu{T, y)) - u{t, 0)5yu(r, y) = 0, 

together with a zero flux boundary condition, dyu{T,0) + u{t,0)'^ = 0. The 
additional left-sided drift contributes to confine the mass in the new frame 
{T,y). It has been proved in [3] that u converges to Ga in large time, where 
Ga is given by 



G„(y) = aexp -ay-^ , / G^{y)dy = u{0,y)dy = M, (12) 



and more precisely, that the following Lyapunov functional L converges to 



In Proposition 11.41 we improve this result with an exponential decay to 
equilibrium. The proof of this result is done as follows. In the new variables, 
the equilibrium state is a gaussian, hence in a linear frame the natural tool 
would be a logarithmic Sobolev inequality established by Gross in [8] that 
we first recall, see [2] for a proof for instance. Although we are dealing 
here with a non linear problem this method will be fruitful. To do so, we 
apply this inequality to a suitable measure, namely Guo{y)dy. It would 
have been natural to apply a logarithmic Sobolev inequality to the measure 








(13) 
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G/3{y)dy but a computation of the entropy dissipation with respect to the 
equilibrium state, G^, leads to the Fisher information expressed with respect 
to Guo{y)dy. Therefore a natural idea was to apply a logarithmic Sobolev 
inequality with respect to Guo{y)dy. 

Lemma 4.1 (Logarithmic Sobolev inequality) . Letv{x)dx = e'yi'p{—V {x))dx 
he a measure with smooth density on Assume that V"{x) > 1 then, for 
u > satisfying u(x)dx = i'{x)dx, we have 



f 

Jo 



2 



On the first hand, let us define the map /3 i— t- C(/3) by 

exp (^-(3y " y ) = M, 

2 

and Vj3 by Va(y) = /3y + ^ + log(7(/3). Such a function satisfies V^'{y) = 1 
and ^ 

C(/3) exp (^-(3y " y) = exp (-^^(y)) = Gp{y), 

hence lemma HTTl can be applied and this yields that 

2H{u\Gf3) < I{u\Gp). (14) 

On the second hand, we recall that the evolution of the relative entropy 
with respect to Gq,, see [3], 

d.,. 1^ . (J(r)-no(r)(l-M))^ 
-L{r) = -I{u\G^,) . 

Furthermore, the relative entropy can be decomposed as follows 

^ f „ lo, (g^) ^ f uir. y) lo. (|M) ) „(, lo. (^ ) 

Recalling the definition of G„g we deduce that 

Guo{y) 



u{t, y) log 



Gaiy) 



dy = J u{T,y)log(-^]dy + J u{T,y)log{exp{{a - UQ{T))y))dy 



Mlog( ^ ) + (q - no(r)) J(r), 



where Gq = C{uo). Therefore, 



2H{u\Gp) = 2Hiu\Guo) + 2Mlog(^] + 2{a - uo{t)) J (t) 
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from which we deduce that 



^L{t) + 2L{t) = 2H{u\Gu,)-I{u\Gu,) + 2M\og(^ 
dr \a 

- 2a{uo{T) - a)(l - M) - {uq{t) - af{l - M). 



Recalhng next a logarithmic Sobolev inequahty (fH|) together with — {uq{t) — 
of {I - M) < 0, it foUows that 

^L(t) + 2L(t) <2Mlog| — ) -2a(uo(T) - a)(l- M). 

CLT V / 

Now, it remains to evaluate the sign of the right-hand side term of the pre- 
vious inequality, which we denote by A. The definition of Co first provides 
that 



Co\ , /Cexp(^-uo(r)y-^jd?/\ / , , G„(y) 



log — = - log ^ = - log / exp {{a - uo{T))y) -^rj—dy 

\^ J~exp(-ay-^j(iy / V./n M 

furthermore since 

-2a(uo(r) - a)(l - M) = 2{a - uo{t)) j ay exp (^-ay - dy 

= 2Mj^"log(exp((a-no(r))y))^|Mdy, 
from Jensen inequality it follows that 

A = 2M (^"log(exp((a - no(r))y))^^dy - log (^"exp((a - no(T))y)^^dy) ) < 
hence 

-^L(r) + 2L(r) < 0. 
dr 

This achieves the proof of proposition 11.41 

□ 

In order to obtain a rate of convergence for the norm we will use the 
Csiszar-Kullback inequality, [S].|12). 

Proposition 4.2 (Csiszar-Kullback inequality). For any non-negative func- 
tions f,g£ such that f^^ f{x)dx = g{x)dx = M , we have that 

11/ - g\\l < 2M r fix) log (l^) dx. (15) 
Jo \9{x)J 
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Corollary 4.3. The following inequalities hold true: 



H{u\Ga) < L(0)exp(-2r), 



u{T,y)-Ga{y)\\i < V2ML(0)exp(-T), 



n{t, x) 



VT+2i 



1 



Remark 5. The 'degradation' of the convergence as M ^ 1 is not contained 
in the rate of convergence but in the intial value L{0). Indeed, the correc- 
tion term to the entropy contains the factor (1 — M) in the denominator. 
Therefore L{0) may become very large as M tends to 1. Hence, even if the 
rate does not depend on the mass, the evolution of L{0) as M — t- 1 affects 
the convergence quality. 

5 Critical mass and trend to equilibrium 

In the critical case, if J(0) < oo, using ([6]), we first notice that the first 
momentum J is conserved. Moreover the stationary states to ([1]) are given 
by the one parameter family: noo(O) exp(— noo(0)a;). Hence Jqo = ^prifi) ^^'^ 
there is only one equilibrium state: fa(a;) = aexp(— ax) with a = 
We also recall the formal computation of the time evolution of the relative 
entropy: 



—H{n\ua) = — n{t,x){dxlogn{t,x) + n{t,0))^dx = —I{n\uo), (16) 



In [3], under the hypothesis that the second momentum is finite, it has been 
proved that the solution to ([T]) converges in relative entropy to fo,. Here, we 
improve the convergence result by precising the speed of convergence. First, 
we give a rate in the case of initially finite second momentum. Then, when 
the third momentum is initially finite, using the HWI inequality we improve 
the rate of convergence. 

5.1 The Lyapunov functional approach 

As we saw in the sub-critical logarithmic Sobolev inequality is a pow- 

erful tool to deal with Gaussian measure. However, in the critical case, the 
stationary state is an exponentially decreasing measure but non-Gaussian. 
In order to use a logarithmic Sobolev inequality, we consider a Lyapunov 
functional approach. Indeed, corrective terms will bring the construction of 
a targeted Gaussian measure. We give a speed of convergence for relative 
entropy with the following theorem. 




where we have used the notation 



fo(i, x)dx = n{t, 0) exp(— n(t, 0)x)dx. 
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Theorem 5.1. Assume that x'^nP{x)dx < +00 and that H{nP\ua) < 
+00, then 

H{n\u^) < -^=L= (^(nVa) + yn°(x)dx) . 

Proof. Let c > be a differentiable function on M_|_. Let us consider the 
nonnegative Lyapunov functional 

F{t) = H{n\ua) + c{t) / —n{t, x)dx. 

Jo 2 

Using equation (fT6]l , we have that 
^F(t) = -I{n\iyo) + c{t)^ (^J^ ^n{t, x)dx^ + c'(t) ^n{t, x)dx. 

We need to bring up a Gaussian measure in Fisher information in order to 
use a logarithmic Sobolev inequality. An easy computation shows that 



d_ 
di 



00 ^2 



n{t, x)dx 



00 ^2 



dx n{t, x)dx lo, 



n(t, x) 
uo{t,x) 



dx 



xn{t, x)dx I loi 



The previous equality (fT7|) leads us to regroup the terms as follows 



(17) 



n{t, x) 
uo{t,x) 



dt 



Fit) 



n{t, X 

) 

n{t, X 



dx log 



00 ^2 



n{t, x) 
i^o{t,x) 

n{t, x)dx 



+ c{t)xdx log 



Mt^x))^"^'^ 4 



dx 



dx log 



n{t, x) 
uo{t,x) 

c{t)x'^^ ' 



+ 2[dx log exp 



' ' ' 5^ log 



i/o(t,a;) 



We define the Gaussian measure GQ^c(t)/2{t,x)dx by 

z/o(t,a;) exp l—£ii}£. 



x)(ix. 



(18) 



G'o,c{t)/2(i,a;)(i3; 



f^uo{t,x) exp 



(19) 



Recalling the definition of the Fisher information together with (|19p , equality 
(fT5|) rewrites as 



n(t, x)(ix. 
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Using next a logarithmic Sobolev inequality (|4.ip for the Gaussian measure 
Go,c(t)/2) it follows that 

Hn\Go^c(t)/2) > c(*)^("|G'o,c(t)/2)> 

hence 

j^F{t) < -c(t)i/(n|Go,.(,)/2) + (c'(t) + ^) y n(t, (20) 

In order to have a rate of convergence, we want to make appear the Lyapunov 
functional on the right-hand side term and we use that 

H{n\GQ^^t)/2) = H{n\va) + I nlog ( — ) 

Jo \Go,c{t)/2j 

which allows rewriting (|20p as 

-F(i)<-c(t)F(t)-c(t) \ n\oA ° +(c^(t)+c(t)^) / —nit,x)dx+c{tf -rn{t,x)a 
dt Jo \Go^c(t)/2/ Jo ^ Jo 4 

(21) 

On the other hand, since J^^ n{t,x)dx = 1, we deduce that 

. f M^) \ , . ( r ( , . c(t)x'^\ , \ 
/ n{t, x) log — - ]dx = log a / exp -n{t, 0)x — dx 

Jo \Go,c(t)/2{t,x)j \ Jq V 4 y y 

f'°° c(t) f°° 

+ (n(t,0)-a)/ xn{t,x)dx + --^ x'^n{t,xlU^] 
Jo 4 Jo 

hence this provides a control on Jq°° n log (^ GQ^"t)/2 } ' ■^^'^^^^ inequality gives 
a first control on the following terms 

. ( r ( r s C{t)x^\ \ , ( r ( , . ^^ c{t)x^\ , , \ 

log I / 6xp ( — n(t, 0)x ^ — I ax I = log ( / exp I (a — n(t, 0))x ^ — I Va\x)dx I 

/■°° c(t) 
> (a — n(t,0)) / xua{x)dx — / x'^Va{x)dx. 







Recalling the conservation of the first momentum, xn{t, x)dx = xvci{x)dx, 
([22]) simplifies as 

' ( \ c(t) f r o , , , r 2 . . . f . ^ , 1 

nlog — > — — / X n[t,x)dx — x Va[x)dx = c[t) / —n{t,x)dx 



Go,c{t)l2j 4 \Jq ' Jq J \Jq A ' 2a2 

coming back to ([2T]). 



-F{t)<-c{t)F{t) + {c'it) + c{tf) / -n(t,x)(irr + |^, 
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hence, choosing c{t) to satisfy 



9 poo o /, \ I ; ' (^2) 

Jg x'^nit, x)ax J 
we have a rate of convergence on the time evolution of F 

jF{t) < -c{t)F{t), (24) 
Using next Cauchy-Schwarz inequahty on the measure n(t, x)dx, we see that 



/ x^n{t,x)da 
Jo 



n{t,x)dx] ( / x'^n{t,x)dx ] > ( / xn{t,x)dx 



hence a first lower bound on c{t) defined by (f23|) : 

1 



c'{t) > -2c{ty, c{t) > 



l + 2t 



For sake of simplicity, we choose c(0) = 1. Using Gronwall lemma in ([2 
we obtain that 

F{t) < F{0) exp 1^ c{s)ds^ < F(0) exp (^-^ log(l + 2t)^ = ■^^=, 

(25) 

Since H(n\i'a) < -?"(*)) tliis achieves the proof. 

□ 
□ 

From theorem 15. II and Csiszar-Kullback inequality ()15p . we immediately 
deduce the convergence 

\\n-iyafi<2H{n\uo,) < ° 



where Cq = V^FiO) = ^2H{'nP\vcc) + x^nO{x)dx. 

Corollary 5.2. Assume that x'^n^{x)dx < +oo and that H{'nP\i'a) < 
+00. For any time t > 0, we have 

II II <r ^0 

\\n — l^aUl < 



(l + 2t)4 

where Cq was previously defined. 

During the proof of Theorem 15.11 we could have found a better lower 



bound for c, defined by (p3|) . and consequently a better rate of convergence. 
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Proposition 5.3. Assume that x'^n^{x)dx < +00 and that H^n^lua) < 
+00. For all f3 < 2/3, there exists a constant C and a time such that the 
following inequality holds true for any time t>tp 

C 

Proof. We have a better lower bound of the second moment with the fol- 
lowing lemma. 

Lemma 5.4. With the previous definition of n and u^, the following in- 
equality holds true 

liminf / x n{t,x)dx > 

Proof. Let {tk)k be the following sequence 

liminf / x'^n{t,x)dx = lim / x'^n{tk,x)dx. 

Recalling Theorem \1.3\ the convergence of n towards Ua holds true. 
Thus, we can find a sub-sequence (ifep)p of {tk)k such that n{tk^., .) converges 
towards almost everywhere. We define 

Up{x) = min(n(tfcp,2;),z^a(x)) < 

Since x'^Up{x)dx < ^ and x'^Up{x) — x'^Va{x) for almost all x, we can 
use Fatou's lemma 

1^00 poo poo poo 

/ x'^Va{x)dx = / liminf x^np(x)(ix < liminf / x^Up{x)dx < liminf / x'^n(t, 

Jo Jo + + Jq t^ + 00 Jq 

□ 
□ 



This lemma provides an upper bound on the right hand side term in (p3 
hence, by definition of the lim sup, for all rj > 3/2, we obtain that 



inf sup 00 2 u ^^ < - 1- 
s>o \ t>s Jq x^n(t, x)dx ) 

Finally, there exists a time t^ such that for t > trj with c(t^) > 0, 

c'{t)>-rjc{tf, c{t) > [ c{s)ds>-log{c{t^r^+r]{t-t^)). 



In the same way as in ()25p . we use the lower bound on c to conclude: for 
any time t > trj 

F{t) < F{t,) -. 

(c(t^)-i+7?(t-t^))^ 
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For all /3 = i < 2/3, the previous inequality proves that (1 + tfF{t) IS 
bounded for t > tp = tr^. 

We conclude this proof by using Csiszar-Kullback inequality (jlSp . 

□ 
□ 



5.2 A better rate 

In the previous paragraph, in order to use a logarithmic Sobolev inequal- 
ity, we have constructed a Gaussian measure with corrective terms. In this 
paragraph we directly use the HWI inequality firstly established in [H] (see 
also [E]), which is adapted for exponentially decreasing measure, and we 
improve the speed of convergence by controlling the second momentum. As 
we have seen in the Lyapunov functional approach, the rate of convergence 
could have been improved by using better estimations on the second momen- 
tum. The Wasserstein distance appears in the HWI inequality and we can 
control this distance with the second momentum. We need a third momen- 
tum uniformly bounded in time for such control on the second momentum. 
This is the following result. 

Theorem 5.5. Assume that sup^g x'^n{t,x)dx < +oo, then, there exist 
to > 0, Ci > and C2 > such that for any time t > to 

1 



Ci+C2t 

Proof. We break the proof into several lemmas. We start by recalling the 
HWI inequality obtained in {14]. see [16] for instance. This inequality binds 
both Fisher information, entropy and Wasserstein distance and can be ap- 
plied to exponentially decreasing measure. 

Lemma 5.6 (HWI inequality). Let v{x)dx = eyi'p{—V {x))dx he a measure 
with smooth density on with V"{x) > and j V {x) e,yi'p{—V [x))dx < 
+00. Then, for n > with finite momentums up to order 2, we have 



n\v'). 

In the particular case of an exponential distribution fo(t, x^dx = n(t, 0) exp(— n(t, 0)x)(ix, 
the HWI inequality reads as 

d Ti \ \ -f^^(^l^o) 

—H{n\ua) = -lin fo) < -tt-j: r- 

dt (n, I/O ) 

This inequality becomes very powerful once it is associated with the following 
inequalities. 
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Lemma 5.7. With the previous definition of n, uq and u^, the following 
inequalities hold true 

1. ff(n|fo) ^ Hlnlua), 

2. W2{n,uof <2j^x^nit,x)dx + ^^. 

Proof. We compare the two terms of the first inequahty 



H{n\vQ)-H{n\va) = I n(t, x) log dx = log ( , ) 0)-a) / xn{t 

Jo \Mx)J \n[t,0)J Jo 

For all X > 0, we recaU that X - 1 > log X then 

H{n\uo)-H{n\ua) = - log ( ^*^^^^ ) + ( ^'^^^^ - 1^ > 0, with ^ xn{t,x)dx = ^. 

Using the definition of the Wasserstein distance ([5]), the second inequality 
is just a consequence of the triangle inequality 



W2in,iyo) < W2{n,5o) + W2{So,i^o) = \l I x^n{t,x)dx+ ^ 



n{t,0)' 



together with (a + 6)2 < 2a'^ + 26^ for a, 6 G R. 

□ 
□ 



The Wasserstein distance is then controlled by the second momentum 

log 2 
Janp 



and n{t,0). Recalling Proposition 13.51 ^or t > = t^, we have a lower 



bound on n{t, 0): 

n(t,0)2 - 

Furthermore using Lemma l5.7t for any time t >tQ, we see that 
—H{n\i'a) < 



dt K + 2 x'^n{t, x)dx ' 

hence after time integration, 

^ > „ , + r r. , nroo\ • (26) 



H{n\va) H{n{tQ,.)\ua) Jt^ K + 2 x'^n{t,x)dx 

Using that X i— )• is convex on M!j_ together with Jensen inequality, for the 
probability measure j^l^f^^f]{s)ds, it follows that 

1 

< 



J^{K + 2 x^n{t,x)dx)j^l[tot]is)ds Jq K + 2 x'^n{t,x)dx t — to 

(27) 



ho,t](.s)ds 
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therefore combining ([26|1 and ([27|) . 



H 



> 



+ 



t-to 



(28) 



Thus, if the second momentum is enough controhed, we obtain a better 
speed of convergence than in the Lyapunov functional approach. 

Remark 6. // x^n{t, x)dx is uniformly bounded in time, then equation 
^M) gives 

H{n\v^) - H{n{to,.)W^)^^^*~*°^- 

However, even if we do not control the second momentum, we can de- 
scribe the behaviour of its Cesaro mean. We can then end with a final 
lemma 

Lemma 5.8. Assume that supj>Q x^n(t, x)dx < +oo, then for any time 
t > to 



t-to 



t fOO 



to JO 



x'^n{s, x)dxds < 



1 



3a(t - to) \ 70 



°° 3 .... 3a 



X n{to, x)dx 



x'^n{to, x)dx I + 



Proof. We start by recalling a simple consequence of Holder inequality 







oo \ 2 

x'^n{t.,x)dx\ < ( / xn{t,x)dx 



x^n{t, x)dx 



with /jjj xn{t, x)dx = ^ then we differentiate and this leads to 



x^n{t, x)dx, 
(29) 



d_ 
di 



x^n{t, x)dx 



x dx(dxn(t,x) + n(t,0)n{t,x))dx = 3n(t,0) / x n(t,x)dx. 

a 



In the same way, we differentiate the second momentum 



d_ 

di 



x'^n{t, x)dx 



x^dx {dxn{t, x) + n{t, 0)n{t, x)) dx = 2—2 



Recalling pip . n(t,0) can be rewritten as 

ad 

"^*'°) = "-2d^ 
allowing us to rewrite equation (fSOj) as 

6 



(30) 

n{t,0) 

a 
(31) 



x'^n{t, x)dx 



d 
'di 



X n{t, x)dx 



a 



3a / X n{t, x)dx+- 



4 dt 



oo \ 2 

x'^n{t, x)dx 
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Finally, we have obtained 

d f f°° , ^ , 3a / /■°° 2 / X , 6 f°° 2 r N , 

— I y X n[t, xjdx — [J X n{t, x)dx I 1= j ^ ' 

hence, after time integration, 

' oo 3a / f°° \^\ /■* f°° 

x^n{t,x)dx — ( / x'^n{t,x)dx\ \ + 3a / x'^n{s,x)dxds 



^ \Jo J j Jto Jo 



x n[tQ,x)dx — / X n[to,x)dx + 



oo 

x^n{to,x)dx — ■ „ ..^_U7~/ 

JO 4 VJo j ) ^ 

therefore, recalling (j29|) . we deduce that 



/ / x^n{s,x)dxds = — I / x'^n(to,x)dx — ( / x^nit^ 

t - to Jto Jo 3a(t / 

□ 
□ 



Consequently, using Lemma 15.81 with (j28p . it follows that 
1 > 1 ^ t - to 



H{n\ua) H{n{to,.)\u^) K + ^ + ^' 



Taking the inverse, this achieves the proof of Theorem 15.5 



□ 
□ 
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